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ABSTRACT 

We present Chandra HETG observations of SS Cygni in quiescence and out- 
burst. The spectra are characterized by He-hke and H-hke K« emission hues 
from O to Fe, as well as L-shell emission lines from Fe. In quiescence, the spectra 
are dominated by the H-like Ka lines, whereas in outburst the He-like lines are 
as intense as the H-like lines. In outburst, the H-like Ka lines from O to Si are 
broad, with widths of 4-14 eV in Gaussian a (1800-2300 km s~^). The large 
line widths, together with line profiles, indicate that the line-emitting plasma is 
associated with the Keplerian disk and still retains the azimuthal bulk motion. 
In quiescence, the emission lines are narrower, with a Gaussian a of 1-3 eV (420- 
620 km s~^). A slightly larger velocity for lighter elements suggests that the lines 
in quiescence are emitted from an ionizing plasma at the entrance of the bound- 
ary layer, where the bulk motion of the optically thick accretion disk is converted 
into heat due to friction. Using the line intensity ratio of He-like and H-like Ka 
lines for each element, we have also investigated the temperature distribution in 
the boundary layer both in quiescence and outburst. The distribution of SS Cyg 
is found to be consistent with other dwarf novae investigated systematically with 
ASCA data. 

Subject headings: accretion, accretion disks — dwarf novae, cataclysmic variables 
— plasmas — stars: individual (SS Cygni) — X-rays: stars 
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Introduction 



A dwarf nova is a close binary system consisting of a weakly magnetized white dwarf 
and a low-mass companion filling its Roche lobe and is characterized by frequent and large- 
amplitude optical outbursts. Mass accretion from the companion occurs via an optically 
thick disk extending down close to the white dwarf, and a so-called boundary layer is formed 
between the white dwarf and the innermost edge of the accretion disk, where the accreting 
matter is braked due to strong friction and finally settles onto the white dwarf. 

After a long debate since the 1970's, it is now widely accepted that the optical out- 
burst arises due to a thermoviscous instability in an outer region of the accretion disk. In 
quiescence, matter from the companion is accumulated in the disk because the hydrogen is 
neutral, and hence the disk viscosity is low. As matter accumulates and the temperature 
exceeds ~10^ K, the viscosity increases due to hydrogen ionization, and matter rushes onto 
the white dwarf at a rate exceeding the supply ra t e frorn the companion. T his latter state is 
the optical outburst (Smalg[l98J; ICannizzolll993l : IOsakilll996l : lLasotall200ll ). In response to 
the phase transition of the outer disk, it is predicted that the physical state of the boundary 
layer changes significantly. In quiescence, the boundary layer is optically thin, because the 
mass accretion rate is low and the radiative cooling is inefficient. The temperature is as 
high as 10^ K, and the boundary layer is a hard X-ray emitter. In outburst, on the other 
hand, the boundary layer can remain optically thick just as for the outer region of the disk, 
because of effic ient radiative cooling. Since its temperature is ^ 10^ K, it emits mainly in 
the E UVhand (jPringlelll977l : iPringle fc Savonijelll979l : lTylendalll98ll : [Patterson fc Raymond 
19851 ). The first observational suggestion to support this phase transition w as obtained from 
a ser ies of European X-ray Observat ory Satellite (EXOSAT ) observations (jjones fc Watson 
1992h and was finally established bv lwheatlev et all J2003h . 



Although the phase transition of the boundary layer is better understood than in the 
past, there still remain several unresolved issues. Although the boundary layer is believed to 
be optically thick in outburst, optically thin hard X-ray emission above 2 keV does not vanish 



completely (INousek et al.l 1 1994 IWheatley et al.l l2003l ) . Unlike the magnetic C Vs in which 



the postshock accretion fiow is essentially one-dimensional, the temperature distribution of 
the boundary layer plasma is not understood very well. The temperature distribution is 
probably tightly linked with the geometry of the boundary layer, for which we have no more 
than a rough sketch. In order to help resolve these issues, we analyzed Chandra HETG 
data of SS Cyg. Since SS Cyg is the brightest dwarf nova in X-rays, it has been observed 
with many X-ray astronomy satellites, such as EXOSAT, Ginga, the Rossi X-ray Timing 
Explorer (R XTE), the Advanced Satellite for Cosmology and Astrophysics (ASCA), Chandra, 



and Suzaku ( lYoshida et al.lll992l : iNousek et al.lll994j : iDone &: Osborndll997l : IWheatley et al. 
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20031 : lMauchJ2004l : Mauche et alJboOSl : llshida et al.ll2007h . Since the Chandra HETG has an 



advantageous energy resolution, we have attempted to approach these unresolved issues by 
means of spectroscopy of He-like and H-like Ka emission lines. We discuss the temperature 
distribution and the geometrical structure of the boundary layer based on our results. In § [2] 
we introduce the Chandra observations of SS Cyg in outburst and quiescence and present 
light curves and spectra. In § [3l the data analysis method and the results are presented. 
The results are discussed in § HJ and the conclusions are summarized in § [51 



2. Observations 



2.1. Chandra data 



The SS Cyg data analyzed in this paper were retrieved from the Chandra Data Archiv^ 
maintained by the Smithonian Astrophysical Observatory. The data in quiescence were 
obtained between 10:28 UT on 2000 August 24 and 00:19 UT on August 25, with a total 
good time interval of 47.3 ks. The outburst data were obtained from 21:09 UT on 2000 
September 14 to 14:15 UT on the following day, with a total good exposure time of 59.5 
ks. Both Chandra ob servations employed the High Energy Transmission Grating (HETG) 
fjCanizares et al.ll2005l ). which allows us to carry out fine spectroscopy with energy resolution 
be tter than CCD reso l ution by rough l y an o rder of magnitude. The se data are analyzed also 
bv lMukai et"al1 J2003h , iMauche et Zl j^OoJ , and lRana et"aD J2006h . By referring to the light 
curve from the AAVSO home pag^, we know that the visual magnitude of SS Cyg during 
the quiescence observation was in the usual range, 11.5-12.5. It then went into a normal 
outburst around September 10. Hence the Chandra outburst observation was performed 
~4-5 days after the onset of the outburst (~3 days after the peak with my — 8.5). The 
visual magnitude during the Chandra observation was my = 9.1-9.7. 



2.2. Light curves in quiescence and outburst 

The HETG consists of two independent grating instruments: the Medium Energy Grat- 
ing (MEG), which covers the 0.4-5.0 keV band with an energy resolution E/AE of ~200 
at 2.0 keV, and the High Energy Grating (HEG), which covers the 0.8-10.0 keV band with 
E/AE ~ 200 at 6.0 keV. In analyzing the data, we used the CIAO software (Version 3.2) 



^http:/ /asc. harvard.edu/cda/ 
^http:/ /www. aavso.org/ 
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and CALDB (Version 3.1.0) provided by the Chandra data center. Figure [T] shows the HETG 
hght curves of SS Cyg in quiescence and in outburst. These hght curves were made by sum- 
ming up the positive and negative first-order photons in the above energy bands with a 512 s 
bin width using the thread dmextract. Backgrounds are not subtracted. The count rate 
of the MEG was 2.90 ± 0.47 c s~^ and 1.00 ± 0.07 c s~^ in quiescence and outburst, respec- 
tively. That of the HEG was, on the other hand, 1.29 ± 0.22 c s"^ and 0.44 ± 0.04 c s"^ in 
quiescence and outburst. The intensity of SS Cyg in outburst is less than that in quiescence 
by a factor of 3 on average and is less variable. The weaker X-ray intensity during optical out- 
burst is associated with the optically thin-to-thick transitio n of the boundary layer when the 
mass transfer rate in the outer disk exceeds a criti cal value (jPringlelll977l : iPringle fc Savonije 



19791 : lTylendalll98ll : [Patterson fc Raymond! 119851). Such even ts are detected by coordinated 



optical, EUV, and X-ray observations (IWheatley et al.ll2003l ). 



2.3. HETG spectra in quiescence and outburst 

In Fig. [2] we show spectra of SS Cyg in quiescence and outburst from the MEG and 
the HEG separately. The spectrum files (PHA files) were extracted from the level 1 data, 
and the positive and negative orders were combined to enhance the signal-to-noise ratio. 
The MEG has a larger effective area in the lower energy range and has clearly detected Ka 
emission lines from O, Ne, Mg, and Si, while the HEG is advantageous at higher energies 
and has detected Ka lines from Fe. Simultaneous detection of the Ka emission lines from 
these elements indicates the existence of optically thin thermal plasma with a wide range of 
temperature distribution, both in quiescence and outburst. In quiescence, the spectra are 
dominated by the H-like lines, whereas in outburst, the He-like lines are as strong as the 
H-like lines. It is remarkable that, as demonstrated in Fig. |2](c) and (d), the emission lines 
in outburst are much br oader than in quiescence. These characteristics are also pointed out 



by iMauche et al.l (120051 ). The profiles of the H-like lines appear rectangular, rather than 



Gaussian-like. 



3. Data analysis 

As described in § 1, we aim to elucidate the nature of the boundary layer of SS Cyg in 
terms of Ka emission lines from abundant metals. We restrict ourselves to those from O up 
to Si, because the quantum efficiency of both the MEG and HEG abruptly drops by a factor 
of ~5 at the iridium M edge ~2.1 keV (Fig. [2]), above which energy the line parameters are 
not very well constrained. Since the MEG has better sensitivity than the HEG in this energy 
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band, we concentrate on the MEG spectra, except for evaluation of the H-hke Si Ka hne, 
for which we also utilize the HEG data (§ 3.3 and § 3.4). 



3.1. Continuum 



In analyzing emission lines, we first attempted to evaluate the continuum spectrum of 
the MEG in the 0.4-2.3 keV band. We masked energy bins, including the emission lines, 
and fitted the remaining spectrum with a bremsstrahlung model undergoing a uniform pho- 
toelectric absorption. Prior to the fitting process, the spectra were rebinned so that each 
energy bin includes at least 100 photons in quiescence and 50 photons in outburst, respec- 
tively. The fits to the quiescence and outburst spectra were marginally acceptable at the 90% 
confidence level (x^ [dof] of 452.5 [441] and 300.7 [241], respectively), with temperatures of 
14.5^3:5 and 4.3lo:? keV and with Nh < 0.3 x 10^° cm'^ and < 1.1 x 10^° cm'^, respectively. 
The fits were improved significantly when we adopted a partially covering absorber with 
(d.o.f.) of 404.3 (440) and 259.1 (240). The temperature of the thermal bremsstrahlung in 
quiescence and outburst is 4:-lto^ keV and I.GIqJ keV, which is covered with an absorber 
with Nh = S.9to,t x 10^^ cm'^ and 4.91JJ x 10^^ cm'^ by 57^^°% and 84+^%, respectively. 
Below we adopt the latter as the continuum model. 



3.2. Emission lines 

Given the continuum model, we retrieved all the energy channels and tried to evaluate 
the line parameters. We fixed all the continuum parameters at the best-fit values in § 13. H 
except for the normalization, and added Gaussians to represent the H-like and He-like Ka 
lines of O, Ne, Mg, and Si. In fitting these Gaussians to the data, all the parameters were 
allowed to vary. Table [H summarizes the best-fit parameters of the Ka lines thus obtained in 
quiescence and outburst, together with the ionization temperature obtained from the He-like 
and H-like line intensity ratio, which are plotted in Fig. [31 

It is clear that the obtained ionization temperatures of all ions in quiescence are system- 
atically higher than those in outburst. They show a wide range depending on the element, 
from 0.3 to 2.0 keV in quiescence and from 0.2 to 1.2 keV in outburst. The emission lines 
seen in the 0.7-0.9 keV bands are from L-shell Fe (mainly Ne-li ke), indicating the existenc e 
of a plasma component with a temperature of ~ 4 x 10^ K (lArnaud fc Raymond! Il992l ). 
These facts unambiguously indicate that the plasma in the boundary layer has a significant 
temperature distribution. The differences in the line intensity ratio distribution, as well as 
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in the line profiles, probably reflect a significant structural difference of the boundary layer 
between quiescence and outburst. 



3.3. H-Iike line profile in outburst 

Unlike the He-like Ka line, which is composed of the resonance, intercombination, and 
forbidden lines, as well as satellite lines, the intrinsic structure of the H-like Ka line is much 
simpler, consisting only of the following two resonance lines: and ^^3/2 

Moreover, their energy separation is only 2.0 eV for Si {E = 2003.9 eV and 2005.9 eV, 
respectively) and is even smaller for the lighter elements. We can thus regard the H-like Ka 
lines of O, Ne, Mg, and Si as consisting of a single component for the HETG resolution, 
and their profile can be used to diagnose the geometry and the dynamical motion of the 
line-emitting plasma. 

From Table [H it is seen that all the H-like Ka lines are significantly broader in outburst 
than in quiescence (see also Fig. [2]). Their widths are in the range 4-14 eV in Gaussian a, and 
corresponding velocities are 1800-2300 km s~^. Their central energies are consistent with 
those in the laboratory. Their profiles, fitted with a broad Gaussian, are shown in Fig. |H We 
then fixed Ec at the laboratory values, which are defined as the emissivity-weighted means 
of the Lya doublet. 

The abscissa axis covers an energy range of 10% of each line's rest-frame energy. The 
best-fit parameters fitted with the broad Gaussian are summarized in table [2l Although all 
the fits are already acceptable at the 90% confidence level, the profile seems to change from 
a Gaussian for lighter elements to a somewhat rectangular structure for heavier elements. 

The rectangular line profile is rem iniscent of a homo geneously emitting thin spherical 



shell that is radially in- or outflowing (IHynes et al.ll2002l ). Of these, an out flow associate d 



with SS Cyg in outburst is detected as the disk wind by the Chandra LETG (jMauchell2004l ). 
The wind is detected as blueshifted absorption lines with a velocity of 2500 km s~^ due to 
intermediately ionized O through Fe, overlaid on a blackbody emission with a temperature 
of 250 ± 50kK. Based on this fact, we infer that the observed H-like emission lines would 
also have been blueshifted, because the plasma below the accretion disk is invisible to us. 
We thus consider the radial outflow picture unlikely. 

Compared with the outflow considered above, a radial inflow is more likely if the in- 
nermost part of the boundary layer becomes optically thin, as in quiescence. The plasma 
in the boundary layer then should take the form of a cooling flow for matter settling down 
onto the white dwarf. The temperature distribution demonstrated in § 13.21 is interpreted as 
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a radial stratification in the cooling flow. The plasma emitting any given line can be treated 
as a geometrically thin spherical shell, since the emissivity of any Ka line is so sensitive to 
the plasma temperature. 

In order to see whether the radial inflow interpretation is acceptable, we have constructed 
a model representing a proflle of an emission line emanating from a geometrically thin uniform 
shell falling onto a white dwarf. The free parameters of the model are the radial velocity v, 
the radius of the shell r^^, the opening polar angle of the shell 60, the inclination angle of 
the orbital plane i, and the central energy of the emission line (see Fig. 

The white dwarf always hides the most blueshifted part of the shell, and the hidden 
fraction is determined solely by Tsh/RwD, where -Rwd is the white dwarf radius. In the case 
that the shell does not expand above the orbital plane enough to form a complete sphere, 
we introduce the polar opening angle 6^ within which the plasma is absent. Some examples 
of the resultant line proflles are shown in Fig. [5]^b). Hereafter we refer to this model as the 
stratifled shell model. 

In fltt i ng the model to the H-like Ka lines, we flxed i at 40°, following the measurement 



of IShafterl (119831 ). We also flxed Ec at the laboratory values. The best-flt parameters of 
the lines with the stratifled shell model are summarized in table [3] and shown in Fig. [61 As 
expected, the fltting to the H-like Si Ka line is improved signiflcantly, because the model 
can represent the shoulders of the proflle appropriately. The fltting to Mg is, however, not 
improved signiflcantly and that to Ne and O became poorer than the Gaussian flt. The 
radial velocity v of all the elements is in the range 3400-4500 km s~^ and does not show any 
sign of slowing down from the Si shell (higher T) to the O shell (lower T). 



Finally, we have tried the diskline model in XSPEC ( lArnaudI Il996l ) to see whether 
the proflle of the H-like Ka lines can be explained by a bulk azimuthal motion of the plasma 
around the white dwarf. The diskline model is intended to represent the emission line 
proflle from a geometrically thin Keplerian disk irradiated by the central ac tive galactic 



nucleus (AGN) via fluorescence, taking the gravitational redshift into account flStellalll990 



Fabian et al.lll989l ). Hence, the parameters are the innermost radius of the disk R^^ in units 
of the Schwarzschild radius and the radial power q of the emissivity e{r) in the form e oc r'^. 
The results of the flts are summarized in table H] and displayed in Fig. [61 The resultant xl 
values of the flts to the O through Si lines with the models introduced so far are shown in 
Fig. [3 

The diskline model can provide a comprehensive understanding to a variety of the 
H-like Ka line proflles from O (Gaussian) to Si (rectangular). On the basis of this result, 
we conclude that the proflles of the H-like Ka emission lines are best interpreted as due to 
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azimuthal bulk motion of the optically thin thermal plasma. 

3.4. Bulk motion of the plasma in outburst 

As demonstrated in § 13.31 the nonrelativistic diskline model is successfully fitted to the 
H-like Ka lines from O to Si in outburst. This implies that the plasma has an azimuthal bulk 
motion. In order to evaluate this bulk motion, we hereafter fit the profiles of the H-like Ka 
lines with two Gaussians that represent the red- and blueshifted line components associated, 
respectively, with the receding and approaching parts of the plasma. The free parameters are 
the average energy of the two line components E^, the energy separation between them Esp, 
the energy widths a, and the normalizations. Among these, E^p is a measure of the absolute 
radial location of the annulus from which the corresponding H-like Ka line emanates, and 
a represents the geometrical width of the annulus. For simplicity, we set the a values, and 
the normalizations of the red- and blueshifted components are constrained to be common. 

The results of the fits are summarized in Table [5] and shown in Fig. [6] and [H with a and 
Esp being normalized by Ec. In addition, we draw confidence contours between Esp and a 
in Fig. [9], since they may couple with each other. 

From Fig.[H]and[ni Esp is found to be smaller for lighter elements, and the 90% confidence 
contours of O and Si are separated, while a does not show any difference among the elements. 
In Table O the rotational velocity v calculated from Esp with v = {Esp/2Ec) c is also listed. 
The rotational velocity of Si is ~1500 km s~^, whereas that of O is only ~10 km s~^. 

4. Discussion 

4.1. Geometry of line-emitting plasma in outburst 

In § l3.3l and § 13.41 we have demonstrated that the profiles of H-like Ka lines from O to Si 
can be explained with the azimuthal bulk motion model (the diskline model), rather than 
the radial inflow model (the stratified shell model) on the basis of fit statistics, and we have 
obtained the bulk motion parameters with the double-Gaussian model. Besides the statistics, 
the high radial velocities ofv = 3400-4500 km s~^ (table [3]) are problematic for the stratified 
shell model. Since the ionization temperature of O in outburst is only ~200 eV, its thermal 
velocity is as small as (^th) = \/SkT/mo ~350 km s~^. Hence the plasma flow is highly 
supersonic, thereby forming a standing shock wave whose temperature kTg = (3/16)/imHW^ 
is as high as 14-24 keV. However, no such high temperature component is detected. The 
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Suzaku observations in outburst, for example, constra in the maximum temperature of the 
plasma to be kT^^^= 8.4 ± 1.0 keV Jlshida et al.ll2007h . 



These results indicate that the bulk motion in the optically thin thermal plasma in 
outburst is azimuthal rather than radial. Although the errors are somewhat large, the 
lighter elements are likely to have smaller Esp (Fig. [8] and [9]). It is natural to expect that 
both the temperature and the azimuthal bulk velocity are gradually reduced as the plasma 
approaches the white dwarf, because it eventually settles down onto the white dwarf. The 
double-Gaussian model is consistent with this expectation in that the lighter elements, which 
have lower ionization temperature, have smaller Egp. 

Since the boundary layer is believed to be optically thick in outburst, it has been unclear 
where the observed optically thin thermal X-ray emission emanates. Our analysis indicates 
that the optically thin thermal plasma also shows a bulk azimuthal motion whose velocity 
decreases as the plasma cools and approaches the white dwarf. We infer that the optically 
thin thermal plasma exists over the o ptically thick disk like an ac cretion disk corona in 
low-mass X-ray binaries, as sketched in [Patterson fc Raymond! (119851 1. 



4.2. Line emission site in quiescence 



In § 12.31 we demonstrated that the emission lines in the quiescent spectra (Fig. [2]) are 
significantly narrower than those in outburst. Nevertheless, the lines still have finite widths 
in the range 1-3 eV in Gaussian a for the H-like Ka lines (Table [1]). They are as large as 
1 part per ~10^ of the line energies E^, and if they are attributed to thermal motion, then 
the following holds: 



a 



1 1 



Er 




m 



where (f^h)^''^ is the root mean square of the thermal velocity. Using the atomic mass number 
A, we obtain 



kT 



mc 



a 



75 I A 

20 



a 



2 eV 



Er 



1 keV 



keV, 



'1^ 



which is 64 keV for O and 51 keV for Si. These temperatures are incompatible with the 
detection of the H-like K a emission lines, b ecause they are emitted from a plasma with a 
temperature of < 1 keV (IMewe et al.lll985l ). Hence, the line broadening should originate 
from a bulk motion whose geometry is so symmetric that the line Doppler shift is observed 
as a line broadening, not as a line shift. O ne may consid er an outflow such as is observed in 
the EUV to soft X-ray band in outburst flMauchdl2004l ). However, the outflow is probably 



- 10 - 



located farther out of the optically thin boundary layer, and most of the receding parts of the 
plasma are hidden by the accretion disk. In f act, the soft X -ray absorption lines originating 



from the outflow in outburst are blueshifted (jMauchell2004j ) 



We consider the observed Ka emission lines from O to Si as emanating from a transition 
region from the optically thick accretion disk to the optically thin boundary layer, where 
the Keplerian motion is gradually converted into thermal motion due to strong friction. The 
average line-of-sight velocity dispersion estimated from the observation is 

which is 620^J?|], 530t\lo, ^^OtHo, and 420±140 km s^^ for O, Ne, Mg, and Si, respectively. 
This is roughly one-tenth of the Keplerian velocity on the surface of the I.2M0 white dwarf 
(^6300 km s~^). Even considering the inclination angle of 40° and the line-of-sight average of 
the azimuthal velocity, (2/7r) sin40°, we obtain the maximum possible line-of-sight velocity 
to be (t')max =2600 km s~^, which can cover the average velocity eq. ([2]) evaluated from the 
observed line width. 

As in outburst, the line width in quiescence can be attributed to azimuthal motion. 
Unlike in outburst, however, where the emission lines can be interpreted as being emitted 
from a cooling plasma, we believe those in quiescence are created in the heating process at the 
entrance to the optically thin boundary layer. In fact, although the statistical significance is 
not very high, the bulk velocity in quiescence estimated from the H-like Ka line widths seems 
to decrease from O to Si. It is likely that the plasma heating results from the conversion of 
the energy of bulk motion. 



4.3. Temperature distribution 

In this section, we attempt to impose a constraint on the temperature distribution of 
the plasma in outburst and quiescence in terms of the intensity ratio of the He-like and 
H-like Ka emission lines. We simply assume that the differential emission measure d{EM) 
is proportional to the power of the plasma temperature. In the following analysis, we adopt 
the spectral model cemekl, in which the temperature distribution is taken into account as 
d{EM) oc (T/Tmax)"c?(logT) oc {T/T^s.x)"'~^dT, where T^ax is the maximum temperature 
of the plasma. Hence, the exponent a obtained b y a fit with this model is larger than the 



true weighting factor by 1 (IDone fc Osbornelll997l ). In the framework of this model, the line 



intensity ratio is uniquely determined once a and Tmax are given. Hence, we can conversely 
give a constraint on these two parameters with the line intensity ratios of O, Ne, Mg, and Si 
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that we have obtained from the HETG data. The allowed parameter regions delineated by 
the 90% errors of the intensity ratios (Table [1]) are shown in Fig. [10] with a different color 
for each element. 

For all the elements, the contours become vertical, and a diverges as T^ax decreases 
and approaches a certain value. This limiting Tmax gives the temperature of the single- 
phase plasma that can reproduce the observed intensity ratio. On the contrary, the contour 
becomes horizontal in the high-Tmax limit. Since the emissivity of atomic lines peaks at a 
finite temperature, the line ratio becomes insensitive to T^ax and is determined solely by a if 
^max becomes too high compared with the emissivity peak temperature. As shown in Fig. [TOl 
there is a parameter region that can be shared by all elements (except for O in outburst). 
Taking into account the constraint Tma x = 8.4 ± 1.0 keV for the outburst obtained from the 
Suzaku observation (jlshida et al.l 120071 ). the common range of a is —0.2 to 0.3 in outbu rst. 



On the other hand, the maximum temperature is 15.4_5 4 keV (IDone fc Osbornd 119971 ) in 



quiescence, which is significantly higher than in outburst, and the common range of a is 
found to be 0.6-1.2 in quiescence. This significant difference is probably the result of the 
geometrical differ ence of the optically th i n bou ndary layer between quiescence and outburst 
(for example, see [Patterson fc Raymond! (119851 )). 



We remark that only the index a of O in outburst is not common with the other 
elements. This is not unreasonable, in view of the shape of the cooling curve of an optically 
thin thermal plasma. Unlike the temperature region T > 2 x 10^ K, the slope of the plasma 
cooling cur ve shows a complicated ch ange as a function of the temperature in the range 
T < 10^ K (IGehrels fc Williamslll993l ). Since a of a given element reflects the slope of the 
cooling curve around the temperature at which the line emissivity of the element shows the 
peak, it is possible in principle that a is different for elements whose emissivity peak appears 
in the range T < 10^ K. 



The index of the temperature distribution a was systematically examined by lBaskill et al 



fl2005[ ) using the complete set of the ASCA nonmagnetic CV data. They found that a for 
the dwarf novae is in the range from 0.7 to 1.8 in quiescence (with a few exceptions) and 
—0.1 to 0.3 in outburst. The a- values obtained fro m our analysis on SSCyg in quiescence 
and outburst are both consistent with their results. Done fc Osbornj (1997) carried out the 
same analysis using the outburst and quiescence data of Ginga and ASCA. Although a in the 
quiescence Ginga data is not constrained very well (2.9+^4), they derived a = OASto.tl 

and 



0.46]'^o;o5 for the normal and anomalous outbursts from the G inga and ASCA dat a, respec- 



tively. The a of the anomalous outburst is outside the range of lBaskill et al.l (120051 ) and that 
of ours. The exponent a may possibly be different between the two types of the outburst 
and may depend on the intensity of the reflection component included in their modeling. 
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The temperature distribution in the postshock accretion column of a magnetic cata- 
clysmic variables (mCVs) can be approximated as T{z) oc z"^^^, where z is the hei ght from the 



white dwarf surface if the cooling mechanism is domi nated by free-free e mission (IFrank et al. 
20021 ). In this case, the cevmekl power a ^ +0.5 (jlshida et al.lll994l ). which lies between 



that in the quiescence and the outburst obtained from this work. 



5. Conclusion 

We have investigated the optically thin boundary layer of SS Cyg both in quiescence 
and in outburst by carrying out high resolution line spectroscopy with the Chandra HETG. 
The spectra contain He-like and H-like Ka emission lines from O to Fe, as well as L-shell 
emission lines from Fe, both in quiescence and outburst. However, the relative intensity of 
the He-like and H-like emission lines of a given element is significantly different; in outburst 
they are nearly equal, whereas in quiescence the H-like lines are much more intense than 
the He-like line. The ionization temperatures evaluated from the line intensity ratio of O, 
Mg, Ne, and Si are in the range 0.3-2.0 keV and 0.2-1.2 keV in quiescence and outburst, 
respectively. 

The other remarkable difference of the line profiles between quiescence and outburst is 
in the line width. The H-like lines from O to Si in outburst are much broader, with widths 
of 4-14 eV in Gaussian a (1800-2300 km s~^), and especially that of Si has a rectangular 
profile, rather than a Gaussian profile. Among a few possibilities, such as a homogeneously 
emitting thin spherical shell that is radially in- or outflowing, an azimuthal bulk motion, 
which can be mimicked by the diskline model, most naturally explains the observed line 
profile. The boundary layer is, however, believed to be optically thick all the way down to 
the white dwarf surface. We thus infer that the line-emitting plasma is located above the 
optically thick disk, like an accretion disk corona commonly seen in low-mass X-ray binaries. 
Since lighter elements tend to have a smaller bulk velocity, the observed line-emitting plasma 
in outburst is a cooling plasma. 

Although less remarkable, the emission lines from O to Si in quiescence also have a finite 
line width with a Gaussian a of 1-3 eV (420-620 km s~^). It is impossible to attribute the 
width to a thermal motion, because the resultant temperature becomes as high as >50 keV. 
Since no red- or blueshift is detected for the line central energies, the radial in- or outflow 
is unlikely. The most natural origin of the bulk velocity is an azimuthal motion. A slightly 
larger velocity for lighter elements, contrary to outburst, suggests that the lines in quiescence 
are emitted in the heating process at the entrance to the boundary layer, where the bulk 
motion of the optically thick accretion disk is converted to thermal energy due to friction. 
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We have also investigated the temperature distribution in the boundary layer within 
a framework of the differential emission measure having a power-law distribution of the 
temperature as (T/Tmax)"- From the line intensity ratio, the index a is constrained in the 
range from —0.2 to 0.3 in outburst and from 0.6 to 1.2 in quiesce nce. These values are 
consistent with those obtained from the previous systematic work of iBaskill et al.l (120051 ). 
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Table 1. Best-fit parameters of a Gaussian to tlie Ka lines and resultant ionization 

temperatures 



E, a I EW /h//hc m 

Ion [keV] [eV] [IQ-^ph. cm-^s"^] [eV] [keV] 



Quiescence 



Si XIII 


1 8649+°-°°^° 


'J--'-'J-3.16 


U.OO_Q 27 


2.9 






Si XIV 


9 nnt^9+o-ooo5 

Z.UUOZ_Q QQ]^2 


9 «n+0-94 


1 c:q+0.24 
i.OO_Q 25 


9.1 


^•'^-1.37 


^•"-'^-0.53 


Mg XI 


1 Qc;9i +0.0024 

1.00Zi_Q QQ20 


A 17+3.83 


n 87+0.26 


2.6 






Mg XII 


-L-4/ZU_Q QQ5g 


9 9c;+1.06 


1 t;4+0.35 
-'-■^*-0.26 


5.2 


1 70+0.90 


-'-•'-'-'--0.15 


Ne IX 


U.yZlU_Q QQQg 


9 ic;+l-21 
^■-■^^-0.95 


9 90+1.03 


3.5 






Ne X 


1 091 O+0.0411 

i.UZiO_Q QQQ3 


-I 01+0.33 


4.68trr^ 


8.7 


2 in+0-76 

^•J-U-1.03 


^•'3^-0.13 


VII 


n c;7n9+0.0399 

U.0(UZ_QQ3gg 


o 9f5+9.34 


Q'yH-5.28 
' -^'-6.34 


6.5 






VIII 


n fic;'?'?+0.0004 

U.UOOO_Q QQQ3 


1 qc:+0.29 
^••^"-'-0.36 


13.78+1^7 


14.9 




'^■"^'-'-0.15 



Outburst 



Si XIII 


1 sfin-^+o-ooii 

i.»bU3_ooo3o 


19 rn+2.49 
1/.0U_^ 00 


4 62+°-^i 
^■'3^-0.39 


46.4 






Si XIV 


Z.UUOO_Q QQ23 


13.88tf0^ 


q 7q+0.34 
"J- '"-'-0.61 


45.4 


0.81+°:}° 


1 99+0.053 
-•^•^^-0.104 


Mg XI 


1 3404+0-0010 

i.04y4i_Q QQ34 


q qt^+1-48 

^•^"-•-i.ee 


t:9+0.57 


25.9 






Mg XII 


1 47^14+0.0014 


10.24+};?^ 


4 09+0-99 

4.Uy_Q 45 


23.4 


0.74+°:- 


74+0.07 
'^-0.04 


Ne IX 


n 01 8+0.0014 
u.yio_o.ooii 


c: 70+1.22 
^- '°-0.97 


1 3 26+2-45 

10.ZD_2 04 


23.2 






Ne X 


1 0194+0-0016 


7 89+1-48 
' •°^-1.34 




25.8 


'-'•^-'--0.22 


49+0-01 

^•^^-0.05 


VII 


n f;7nQ+o-ooi7 


4 42+0-96 

^•^^-0.81 




97.5 






VIII 


n fi5sn+o-ooo7 

U.DOOU_Q QQQg 


qn+0.68 


47 -, t:+9.23 
^'•-■-^-5.70 


65.7 


54+0-17 
'-'•^^-0.09 


n 99+0.02 

^•^^-0.02 



Table 2. Best-fit parameters of a Gaussian to the outburst H-like Ka lines 



Ion 


[keV] 


cr 

[eV] 


/ 

[10-^ cm-2s-i] 


a/E, 
[10-3] 


xVd.o.f. 


Si XIV 
Mg XII 
Ne X 

VIII 


2.0052 
1.4724 
1.0218 
0.6535 


7 c:9+1.69 
' -0^-1.50 

3.481^7^ 


74+0.41 
■J- ' ^-0.39 

4.39t°:™ 

11.75l}9l 

48.371?:°} 


7 1 c;+0.74 
' •-'-O-0.82 

7.561}:^° 

7 qc:+1.65 
' ■'-^'-^ — 1.47 

cr qq+1.62 


31.6/35 
19.7/29 
21.0/18 
6.1/7 



Fixed at the laboratory values. 
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Table 3. Best-fit parameters of tlie stratified shell model to the outburst H-like Ka lines 



Ion [keV] [deg] [km s^"^] 



Si XIV 


2.0052 


651^2 


9 1 +00 


39201^^0 


15.5/33 


Mg XII 


1.4724 






40301^7° 


19.0/27 


Ne X 


1.0218 


441}^ 


1 q+0.4 


45101^^^° 


16.9/16 


VIII 


0.6535 




1 q+0.6 
-'-■'-'-0.2 


338011° 


6.1/5 



Note. — Inclination is assumed to be 40°. 

Fixed at the laboratory values. 

^ Opening polar angle of the shell within which the plasma 
is absent. 

Ratio of the radius of the shell to that of white dwarf. 
Radial velocity of the shell. 



Table 4. Best-fit parameters of the diskline model to the outburst H-like Ka lines 



Ion 


E,- 
[keV] 




-Rin 

[10^ cm] 


/ 

[10--^ cm-2s-i] 


xVd.o.f. 


Si XIV 
Mg XII 
Ne X 

VIII 


2.0052 
1.4724 
1.0218 
0.6535 


q qi -1-0.44 
'^•'^-'--0.74 
9 yq+O.SG 
^- ' 'J-0.55 
9 -1 cr-fO.lS 
^•-■-^-0.23 

i.JO_g 23 


' •U-2.4 
9 .+1.4 
^■^-0.9 

2.5_2 2 


3.73 ±0.39 
4 26+°-^^ 

X2.55 2 13 

53.601^;^^ 


15.4/34 
17.4/28 
19.6/17 
2.3/6 



^ Fixed at the laboratory values. 

^ Radial power of the line emissivity e oc r'^. 
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Table 5. Best-fit parameters of tlie double Gaussian model to tfie outburst H-like Ka lines 







/ Ec Esp" 1 Ec 


V sin i'^ 


xVd.o.f. 


Ion 


[keV] 


[xlO-3] [xlO-3] 


[km / sec] 





Si XIV 


2.0052 




10.6 ± 1.1 


1590 ± 170 


19.4/34 


Mg XII 


1.4724 


4.6l?:« 


q q+2-3 


1490^^^° 


18.0/28 


Ne X 


1.0218 




< 8.6 


< 1290 


21.1/17 


VIII 


0.6535 




< 7.3 


< 1100 


6.1/6 



^ Fixed at the laboratory values. 

^ Common width of the blue- and redshifted components in 
Gaussian a. 

^ Energy separation of the blue- and redshifted components. 

Average line-of-sight velocity of the blue- or redshifted line 
component calculated from Esp/E^. 
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Fig. 1. — HETG light curves of SS Cyg in quiescence and outburst with a bin size of 512 s: 
(left) MEG light curves, (right) HEG light curves. The origin of the time axis coincides with 
the observation start time both for quiescence and outburst. 
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Fig. 2. — Spectra of SS Cyg in quiescence and outburst from (a) the MEG and (b) the HEG. 
Blow-up of (c) the MEG and (d) HEG spectra around the He-hke and H-hke Si Ka hues. 
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Fig. 3. — The ionization temperatures obtained from the ratio of He-hke and H-hke Ka hues 
for each element. 
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Fig. 4. — Results of a Gaussian fit to the MEG spectra of the H-hke Ka hues. The HEG 
spectra are also used for the Si Ka line. The fits are acceptable at the 90% confidence level. 




E / E„ 



Fig. 5. — (a) Schematic view of the stratified shell model, and (b) an example of the line 
profiles. The adopted parameters are i =40°, rsh/-RwD = 1-45, and v = 3000 km s~^, and 
the line profiles are shown at opening angles 6o = 0°, 25°, 50°, 75°, and 85°. 
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Fig. 6. — Results of the fits of (top), the stratified shell model (middle) the diskline model, 
and (bottom) the double- Gaussian model to the H-like Ka lines. 
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Fig. 7. — xl values of the fits to the O through Si fines witfi tfie single- Gaussian model, tfie 
stratified shell model, and the diskline mode. 
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Fig. 8. — (a) Resultant width a and (b) the separation energy E^p of the two Gaussian line 
components in the double-Gaussian model, both normalized by the line central energy Ec- 
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Fig. 9. — The 90% and 99% confidence contours (dotted and continuous lines, respectively) 
between the line separation energy Esp and the line width a in the double-Gaussian model 
for the H-like Ka lines from O to Si. Esp decreases from Si to O, and the 90% confidence 
contour of Si is isolated from that of O. 




[keV] 



Fig. 10. — Areas allowed by the observed line intensity ratio of O, Ne, Mg, and Si in the 
plane of the two cemekl parameters a and T^ax- The allowed areas are defined with the 90% 
confidence contours. Left and right panels are those for quiescence and outburst, respectively. 
Vertical dotted hues at 15.4l|j_and 8.4 ± 1.0 keV are the values of T^ax from llshida et al.l 



( 120071 ) and lDone fc Osbornd (Il997l ). respectively. 



